Impurity States in Graphene 
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Defects in graphene are of crucial importance for its electronic and magnetic properties. Here 
impurity effects on the electronic structure of surrounding carbon atoms are considered and the 
distribution of the local densities of states (LDOS) is calculated. As the full range from near 
' field to the asymptotic regime is covered, our results are directly accessible by scanning tunnelling 

' microscopy (STM). We also include exchange scattering at magnetic impurities and eludicate how 

strongly spin polarized impurity states arise. 
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Graphene, a recently discovered allotrope of carbon 
and the first known example of a truly two-dimensional 
(2D) crystal ^ Ja] has unique electronic properties 
H H H H 11 li [H m m El El such as an 
exotic quantum Hall effect with half-integer quantization 
of the Hall conductivity |3, , finite conductivity at zero 
charge-carrier concentration [^, strong suppression of 
weak localization E3 1 ^^c. The peculiar 2D band struc- 
ture of graphene resembles ultrarelativistic electron dy- 
namics near two nodal points in the Brillouin zone. This 
provides a new bridge between condensed matter theory 
and quantum electrodynamics (index theorem and the 
half-inte ger quantum Hall effect Q , relativistic Zitterbe- 
wegung Il6| and the minimal conductivity E3l> "Klein 
paradox" |l7| and anomalous tunnelling of electrons in 
graphene through potential barriers ilS]). Unexpectedly 
high electron mobility in graphene and its perfect suit- 
ability for planar technology makes it a perspective ma- 
terial for a next-generation, carbon-based electronics ^1] . 

Impurity states are important contributors to these un- 
usual properties. Graphene is conducting due to carriers 
that can either be introduced by a gate voltage P, H 
or by doping P, H 0|. This situation is very reminis- 
cent of doped semiconductors, where the desired prop- 
erties are obtained by creating an impurity band. Re- 
cent progress in scanning tunneling microscopy (STM) 
made it possible to image impurity states for a wide class 
of materials with very high spatial resolution. This so- 
called "wave function imaging" yields local images of the 
impurity-induced wave function. Examples of wave func- 
tion imaging range from unconventional superconductors 
0,01 to semiconductor s l2^l2ll . magnetic metals ^ 
and graphite surfaces 0, 124l25j . It allows one to inves- 
tigate the formation of the impurity band and the asso- 
ciated electronic properties. Theoretical modelling and 
STM measurements of near impurity site effects can be 
compared and thus eludicate, e.g., magnetic interaction 
mechanisms [21I Eij . 

The purpose of this letter is to address the question 
of electronic properties of single and double impurities 



in graphene in connection with future STM experiments 
and impurity-induced ferromagnetism. The impurity 
states are characterized by their energy and by their 
real space wave functions that determine the shape of 
the resonance. In contrast to previous studies |3, we 
consider the real space structure of the electronic state 
in the range from the impurity site to the asymptotic 
regime, its dependence on the potential strength and the 
spin exchange interaction. We will address the impu- 
rity states as an input into a subsequent discussion of 
impurity-induced bands. 

The honeycomb arrangement of carbon atoms in 
graphene can be described by a hexagonal lattice with 
two sublattices A and B. (See, e.g., jl^). With the 
Fermi operators Ci and di of electrons in cell i at sub- 
lattice A and B, respectively, we describe a single and 
two neighbouring impurities by V"s = J/ocJcq and Vd = 
t^o(cJco-|-dJdo) + f^i(cJrfo+rfoCo)- Here Uq is the potential 
strength and Ui the change of sublattice hopping between 
the two impurity sites. Related to the current research 
are questions about imp urities in graphite that have been 
studied with STM [2g. Only the atoms above hollow 
sites are seen in STM on graphite. We find that im- 
purity states in graphene are qualitatively different from 
those in graphite because of the sublattice degeneracy 
that is reflected in a complicated sublattice structure of 
impurity induced resonances. 

We find that impurity scattering produces low energy 
resonances with the real space structure and the reso- 
nant energy i^imp as function of Uq (and Ui) clearly dis- 
tinguishing between single and double impurities. For 
single impurities we find in agreement with Loktev |lO|. 
that -Eimp is well described by 
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where W is the band width. Hence the resonance energy 
E\ai-p approaches zero for Uq 00. Only strong single 
impurities (i.e. C/q 10 eV) are capable of producing 
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resonances within leV of the Dirac point. This result is 
similar to the impurity states observed in unconventional 
superconductors with Dirac spectrum [isj . 

The resonance of a double impurity is basically deter- 
mined hy Uq — Ui. Its energy coincides with the Dirac 
point at finite Uq — Ui = 3t, where t « 2.7 eV is the 
nearest neighbour hopping parameter of graphene. 

We give a detailed description of the local density of 
states (LDOS), the real space fingerprint of impurities in 
graphene, both near the impurity and in the asymptotic 
regime at large distances. Near the impurity site that 
LDOS exhibits an intricate pattern. A single strong im- 
purity placed on one sublattice produces a peak in LDOS 
at low energies that is large on the other sublattice. At 
large distances these impurity resonances have wave func- 
tions V' that asymptotically decay as oc 1/r. 

We will consider potential scattering (nonmagnetic) as 
well as a magnetic impurities, i.e. spin dependent scat- 
tering. In the latter case the impurity induced resonance 
will exhibit a spin dependent splitting that might lead 
to a strong spin polarization of the impurity state. This 
observation, we believe, is important for the the discus- 
sion of moment formation and possible magnetic order in 
graphene. 

To start with our theoretical models we describe the car- 
bon pz-electrons within the tight binding approximation 

- c{k) - 
d{k) 

r(fc) ^^^^ ^ tE-^ie^'^^"^-^^'. 

r^B denotes the Brillouin zone volume and c(fc) (c?(fc)) 
are the /c-space counterparts of Ci (d.i).^l\ 

The full Green's function in real space G{i,j,E) will 
be obtained using the T matrix formalism: 

G(z, J, E) = GO(z - 3, E) + G\i, E)T{E)G'>{-j, E). (2) 

Therefore the unperturbed Green's function G'^{i,E) in 
real space is calculated from its k-space counterpart 
(^{k, E) = [E — Hk + iSe)^"^ by Fourier transformation. 
Numerical problems in carrying out the Fourier integrals 
are avoided, by linearizing the bandstructure in a vicinity 
of the Dirac points, where all singularities occur. Out- 
side these regions the full tight-binding bandstructure is 
taken into account. Finally the T matrix is given by 



by = L 



with -^{k) = 

where S,{k) 



and 



T(^)= (1- 14(d) GO (0,^) 
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with v; = [/o ^ Q Q j and Vd = ^ [/g j being the 

impurity potentials in k-spacc and matrix form. Poles of 
the T-matrix corresponding to impurity resonances oc- 
cur, if det(l - ys,dG°(0, E)) = 0, i.e 

C/oG°i(0,^)-l = (4) 

for a single scatterer and 

(1 - C/oG?! (0, £;))2 - U^Gl, (0, i?)G?2 (0, £;) = (5) 




Figure 1: (Color online) The energy -Bimp of the impurity reso- 
nance as function of the potential strength Uo is shown for sin- 
gle impurities and double impurities with f/i = 0. For the sin- 
gle scatterer iSimp obtained from our tight binding calculation 
is compared to the result obtained from the fully linearized 
bandstructure (Eqn. with fitted bandwidth W = 6.06 eV. 



for double impurity with Ui — Q — this case, we re- 
fer to as scalar impurity. Near the Dirac points we 
have |Re(G°)| 3> |Im(G°)|, so that the impurity reso- 
nances i?inip as function of Uq can be calculated from 
the previous two equations considering only the Re(G°): 
The resulting real impurity energies as a function of 
Uq are shown in figure ^ Adjusting the bandwidth 
parameter W in Ecyr. ^ to fit our £'imp(C^) yields 
W = 6.06 ± 0.02 eV. I.e. W differs shghtly from the 

bandwidth W = ^vi^J^ « 6.3 eV obtained by assum- 
ing linear dispersion in the entire Brillouin zone, [loj 

It is quite remarkable that a pair of neighbouring 
scatterers produces a resonance at the Dirac point for 
Uq — 3t fa 8.1 eV, while for a single impurity this occurs 
only in the limit of infinite potential strength. This effect 
can be attributed to the existence of two nonequivalent 
Dirac points in the Brillouin zone. As a consequence at 
E = Q the onsite Green's function G'^{0,0) has finite off 
diagonal G?2(0,0) = G^i(0,0) = but vanishing di- 
agonal components resulting via equations^and|Slin the 
characteristic i?imp(^^o) curves. 

For double impurities with sublattice hopping change 
Ui, it follows directly from the secular equation, that the 
impurity energy as a function of Uq and [/i is obtained 
from the scalar case by replacing Uq with Uq — Ui. 

We obtained the LDOS N{r,E) 
— ifm 1^^. ^. $i(r)G(i, j, i?)<i>j(r)^ in presence of 
impurities as a function of position and energy, where 
<i>j(r) = [(j)'i{r),(j)f{r)] with (j)l"'\r) being carbon p^- 
orbitals located in the unit cell i at sublattice A and 
B respectively. This LDOS of impurity resonances at 
^-imp = — O.leV for a single and scalar impurities are 
shown in figure |21 (22j The formation of virtual bound 
states (VBS) due to impurity scattering is clearly visible. 

These VBS are a general feature of localized states hy- 
bridizing with a continuum of delocalized states. They 
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Figure 2: (Color online) r-dependent LDOS &t E = -Bimp = 
— O.leV for a single impurity with Uo — 45 eV (left) and for 
a scalar double impurity with Uo — 6.9 eV (right) encoded 
corresponding to color bar. The impurity sites are marked as 
big red dots in the center of the images. 



have been observed in many systems ranging from ele- 
mentary metals [s^l to d-wave superconductors [2^ . De- 
tails of the real space image are however system specific. 
Here, the threefold (-D3/1) symmetry of the VBS for a 
single impurity and twofold {D2h) symmetry for a dou- 
ble impurity are direct consequences of the lattice sym- 
metry. The {D^h) symmetric single impurity state re- 
sults in sixfold symmetric Fourier transformed scanning 
tunnelling spectra 31]. Furthermore the peculiarities of 
the bandstructure of graphene manifest themselves in the 
near field characteristics of impurities: A single impurity 
in sublattice A induces an impurity state mostly local- 
ized in sublattice B and vice versa due to the fact that 
G\^{i,E) < G\2{hE),GlS,E) for E which can 

be attributed to the existence of two nonequivalent Dirac 
points as explained above. 

The site projected DOS N{i, E) can be obtained from 
the full Green's function N{i,E) = ~l/TTlmG{i,i,E), 
where each of the diagonal matrix elements corresponds 
to one sublattice. For the single impurity and Uq from 
10 eV to 40 eV the LDOS at the impurity, nearest neigh- 
bour (NN) and next NN sites are shown in figure 01 One 
sees that for vacancies with lOeV^ ^ Uq 20 eV p, but 
not for weaker potentials, an impurity resonance should 
be clearly observable within 1 eV around the Dirac point. 

It further illustrates the localization of the impurity 
state on sublattice B, when the impurity is in sublat- 
tice A, as well as the reduction of LDOS at the impurity 
site for strong repulsive potential. The double impurity 
respects pseudospin symmetry and is much more sensi- 
tive to weaker potentials as can be seen from figure ^ 
Clearly, Ut — Uq ~ Ui is the most important parameter 
determining the shape of LDOS in the case of a double 
impurity: The results for {Uq = 4eV, Ui = — 2eV) and 
{Uq = 6eV, Ui = OeV) are virtually indistinguishable at 
the impurity resonance but differ slightly below it. For 
large distances r 3> hv[/Einip from the impurity site we 
obtain for the changes in LDOS AN{r, Ei„^p) cx 1/r in 
agreement with ^3,^3 for EtU considered types of impu- 
rities. Note the contrast to a single hard-wall impuritiy. 
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Figure 3: (Color online) LDOS at the impurity site (left), NN 
(middle) and NNN sites (right) is shown for a single impurity 
with potential Uo = 10 eV (a), 20 eV (b) and 40 eV (c). 



0.25 
0.2 

0.15 
0.1 

0.05 








0.12 
0.1 

0.08 
0.06 
0.04 
0.02 






0.14 
0.12 
0.1 






(a) 
(b) 
(c) 


i 1 

\l \ 

^ Al 




V -f i 




0.08 
0.06 
0.04 
0.02 



























-1 -0.5 0.5 



1 -1 -0.5 0.5 1 
Energy relative Fermi level [eVJ 



1 -0.5 0.5 1 



Figure 4: (Color online) As figure |3 but for scalar double 
impurity with Uo — 4eV (a) and 6eV (c) as well as for double 
impurity with Uo ~ 4eV and Ui = — 2eV (b). 



i.e. Uo — > 00, with 1/r^ asymptotics of AN{r, Eiaip) 

If the impurities have a magnetic moment, exchange 
scattering of the graphene pz-electrons and the spin S 
localized at the impurity site will occur. As long as 
the exchange coupling J does not exceed a critical value, 
Kondo screening of the spin S by the band electrons can 
be neglected and the impurity spin acts as local magnetic 
field: The effective scattering potential is rcnormalized 
to Uq ± J. The resulting change in spin-polarized (SP) 
LDOS in the vicinity of a single impurity is shown in 
figure |S1 for Uq — 12eV and J = 2eV. The exchange 
splitting of the resonances in the two spin channels is 
approximately O.lSeV. This type of exchange scattering 
also affects decay lengths and oscillation periods of the 
induced spin density variations and therefore provides a 
possible mechanism for long range exchange interactions. 

For double impurities the effect of exchange splitting 
is much more pronounced within a realistic parameter 
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Figure 5: (Color online) SP-LDOS at the impurity site (left) 
and a NN site (right) is shown for a single magnetic impurity 
with Uo = 12 eV and J = 2 eV. 
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Figure 6: (Color online) SP-LDOS at NN site of a double 
impurity with Uo = 5eV, Ui = -2eV and J = leV (left) 
and 2 eV (right) are shown. 

range: As figure shows, exchange scattering can pro- 
duce strongly spin polarized impurity states. The impu- 
rity resonances of one spin channel can be pushed close 
to the Dirac point or the impurity levels are split even 
below and above it. Depending on the type of impurities, 
the spin polarization of the impurity states can strongly 
depend on doping: In the example with J = 2 eV, the 
VBS above the Dirac point can be occupied by spin down 
electrons due to n-doping. 

It was demonstrated recently js^l that ferromagnetism 
of sp electrons in narrow impurity bands can be charac- 
terized by much higher Curie temperatures than those 
typical for traditional dilute magnetic semiconductors. 
The reasons are a suppression of T-matrix renormaliza- 
tion of the Stoner exchange parameter and a high magnon 
stiffness constant. So the Curie temperature is basi- 
cally determined by single-particle Stoner excitations, in 
a sharp contrast with d-electron magnets. Hence the 
impurity band associated with the magnetic impurities 
considered in this letter can be a promising candidate 
for facilitating high temperature ferromagnetic order in 
graphene. 

In conclusion, we have calculated the LDOS of impu- 
rity resonances in graphene from the near field to the 
regime of asymptotic 1/r decay. The near field LDOS 
are directly observable by STM and comparison of up- 
coming experiments with our predictions will eludicate 
the nature of impurities in graphene. We also find that 
impurity resonances in graphene are very different from 
the impurity states observed in graphite because of the 
two sublattice structure in graphene. 

We showed further how spin polarized impurity states 
can result from exchange scattering at magnetic impu- 
rities and their sensitivity to doping. The resulting for- 
mation of spin polarized impurity bands may give rise 
to long range exchange interactions and magnetic order, 
that can be directly studied by spin-polarized STM. 
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